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ABSTRACT 

Aims. We investigate the suitability of several molecular tracers, in particular CN N=2-l line emission, as a confusion-free probe of 
04 ' the kinematics of circumstellar disks. 

Methods. We use the IRAM 30-m telescope to perform a sensitive search for CN N=2-l in 42 T Tauri or Herbig Ae systems located 
principally in the Taurus-Auriga region. 13 CO J=2-l is observed simultaneously to provide an indication of the level of confusion 
with the surrounding molecular cloud. The bandpass also contains two transitions of ortho-H 2 CO, one of SO and the C 17 J=2-l line 
T ^ ' which provide complementary information on the nature of the emission. 

Results. While 13 CO is in general dominated by residual emission from the cloud, CN exhibits a high disk detection rate > 50 % 
in our sample. We even report CN detection in stars for which interferometric searches failed to detect 12 CO, presumably because 
O [ of obscuration by a foreground, optically thick, cloud. Comparison between CN and o-H 2 CO or SO line profiles and intensities 

divide the sample in two main categories. Sources with SO emission are bright and have strong H 2 CO emission, leading in general 
to [H 2 CO/CN]> 0.5. Furthermore, their line profiles, combined with a priori information on the objects, suggest that the emission is 
coming from outflows or envelopes rather than from a circumstellar disk. On the other hand, most sources have [H 2 CO/CN]< 0.3, no 
SO emission, and some of them exhibit clear double-peaked profiles characteristics of rotating disks. In this second category, CN is 
likely tracing the proto-planetary disks. From the line flux and opacity derived from the hyperfine ratios, we constrain the outer radii 
of the disks, which range from 300 to 600 AU. The overall gas disk detection rate (including all molecular tracers) is ~ 68%, and 

■ decreases for fainter continuum sources. 

Conclusions. The current study demonstrates that gas disks, like dust disks, are ubiquitous around young PMS stars in regions of 
1 isolated star formation, and that a large fraction of them have R out > 300 AU. It also shows the potential of the CN N=2-l transition 

, to probe the kinematics of these disks, thereby measuring the stellar mass, using high resolution observations with ALMA. 

^"T* , Key words. Stars: circumstellar matter - planetary systems: protoplanetary disks - Radio-lines: stars 
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■ 1 . Introduction stars with disks are generally located in molecular clouds whose 

optically thick 12 CO emission may mask the disk emission of the 

; Masses of young stars are a fundamental property with impacts star xhe second limitation was technical: it was difficult to find 

> . on theories of the stars' formation, calculations of their evo- young stars with masses smaller than 5 Mq that have detectable 
lution to the main sequence, and the chronology of the plan- gaseous circumstellar disks. A possible solution to the first prob- 

. ets that may be, or will be, associated with them. Dynamical lem could be obser vation in the lines of 13 CO or C ls O because 
g . techniques provide the only absolutely reliable measurements the optical depth of the moleC ular cloud would be less in these 
of mass. Because most binary orbits are at least a few years transitions than in i2 CO However, the disk emission in the lines 
long, the quickest way to accomplish a dynamical measurement of the CQ isot0 p i ogues would of course be less than that of the 
is by mapping the Keplerian rotation of the circumstellar disk i2 CO lines . We judged that, at the time of our observations, the 
of a young star ( |Gmlloteau & Dutrey| | 1998t | Simon etalj | 2000|) . instrume ntation available to us had insufficient sensitivity to ob- 
These early stodies^and especially the more recent observations tain mass measurements to the required precision, < 5%. The 
by I Schaeferet aU d2Q09D, were all carried out by mm-wave in- commissioning of the Atacama Large Millimeter/Submillimeter 
terferometry in the J=l-0 and J=2-l transitions of 12 CO and re- AlTay (ALMA) offers the promise to remove both obstacles to 
vealed two limitations of the method. The first is intrinsic: Young mass measU rement of young low mass stars by mapping the ro- 
tation of their disks. 
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Send offprint requests to: S. Guilloteau, e-mail: 

guilloteau@obs.u-bordeauxl.fr Our goal in this work is to assess the applicability of sev- 

* Based on observations carried out with the IRAM 30-m telescope, eral molecular tracers for confusion-free and high sensitivity 

IRAM is supported by INSU/CNRS (France), MPG (Germany) and measurements of disk kinematics using ALMA. After 12 CO, the 

IGN (Spain). strongest lines in disks are the rotational transitions of 13 CO J=2- 
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1 or 3-2, HCO + J=3-2 and the hyperfine split lines of CN N=2- 
1 (iDutrev et alJll997h lOberg et alJl2010h . Howev er, 13 CO line 
emiss ion can also be affected by confusion (e.g. IDutrev et alj 
I2003L for BP Tau). The situation is less clear for HCO + J=3- 
2, but this transition can still have substantial optical depth and 
be ex cited in clouds (see for example CW Tau in Salt er et alj 
1201 ll) : furthermore, its relatively high frequency requires bet- 
ter observing conditions. Other previous searches for molecular 
tracers in disks were conducted with the single-dishes (princi- 
pally the IRAM 30-m and JCMT) or mm/sub-mm arrays (IRAM 
PdBI, SMA and CARMA). The former were limited in sensitiv- 
ity (e.g. lSalteretal]|20Tll for a study of CN N=2-l and HCO + 
J=3-2), while, because of the need for suffi cient UV coverage , 
the latter could only observe a few stars (e.g. lOberg et al.l l2010). 

Recent advances in receiver and back-end technology at the 
IRAM-30m telescope encouraged us to carry out a new assess- 
ment of molecular line tracers in circumstellar disks, particularly 
the CN N=2-l line hyperfine split lines. We report here our sys- 
tematic search for CN N=2-l emission in 40 young stars located 
in the Taurus/Auriga star forming region (SFR). Together with 
previous observations, 46 stars in this SFR have now been ob- 
served with high sensitivity in this transition, a gain of a factor 
almost four in sample size. §f2] describes our sample, the obser- 
vations, and their analysis. §[3] presents our results for the CO 
isotopologues, CN, H2CO, and SO molecules. §j4]discusses our 
assessment of CN as a disk tracer and comments on stars with 
outflows and on general properties of the observed line emission. 
$5] summarizes our results. An Appendix shows all the observed 
spectra and provides details for many of the stars. 



2. Observations and Data Analysis 

2. 1 . Source Sample 

We focus here principally on the Taurus Aurigae region, to which 
we add CB26, an embedded source in an isolated Bok glob- 
ule, and HD 163296, an isolated HAe star. Our study is thus 
only relevant to regions of isolated star formation. Table Q] lists 
the stars involved in this study; the 40 stars observed anew at 
the IRAM 30-m radiotelescope appear without specific sym- 
bols. This sub-sample was observed in an homogeneous way. 
It is completed by observations on a few sources (marked by 
specific symbols in the Table) obtained independently, also with 
the IRAM 30-m. 13 CO, CN and the ortho-H 2 CO 3 B - 2 12 tran- 
sitions were observed towards the Herbig Ae star HD 163296, 
together with the S O 65 - 54 line. Data from DM Tau is from 
Dutr ev et alJ (Il997l) . The intensity reported for LkCa 15 is from 
a noisy IRAM 30- m spectrum, bu t is con sistent with the inter- 
ferometric data of Chap illon et al.l d2012h . CI Tau, CY Tau and 
GO Tau were detected with the IRAM 30-m, and late r imaged 
at the IRAM Plateau de Bure array dGuilloteau et alj|2013l in 
prep). The overall sample cov ers a wide range of sp ectral types 
(from M5.5 to AO, taken from Luhman et alI (l201C )), luminosi- 
ties (0.2 to 45 L Q ), and 1.3 continuum flux densities (from < 4 
to ~ 800 mJy), but is not complete with respect to these quan- 
tities. No specific account about extinction and location with 
respect to molecular clouds was made when selecting sources. 
Sources with outflows were deliberately included. The epoch 
2000 coordinates are from interferometric mm continuum im- 
ages dGuilloteau et alj|201 lb when available, or SIMBAD other- 
wise. 



2.2. Observations 

Most observations were carried out with the IRAM 30-m tele- 
scope from Dec 4 to Dec 8 201 1, under excellent weather con- 
ditions (precipitable water vapor content below 1 mm). A few 
additional sources were observed on the afternoon of May 1st, 
2012 using the same instrumental setup. The weather was good 
(water vapor content about 3 mm), but a few clouds occasionally 
produced anomalous refraction. The absolute calibration of this 
second set is thus more uncertain. 

The dual-polarization, sideband separating receivers were 
tuned to cover approximately 7.8 GHz of bandwidth around 
224.3 GHz in the upper sideband (USB) and 208.7 GHz in 
the lower sideband (LSB). The 8 units of the FTS back-end 
(each of them being made of 3 adjacent sub-units) was cov- 
ered more precisely 4 frequency ranges (with 2 polarization 
each) at [220.325,224.365] and [224.055,228.105] GHz in USB, 
and [204.655,208.695] and [208.375,212.425] GHz in LSB. The 
back-end provides a spectral resolution of 195 kHz, which trans- 
lates into 0.25 - 0.28 km.s~' in this frequency range. The spectral 
setup includes several lines (listed in Table |2| which can be rel- 
evant for disks, in particular all 19 hyperfine components of CN 
N=2- 1 . The sideband separating mixers provide a mean rejection 
around 13 dB. However, the presence of strong lines spread over 
the band indicate some variation across the band. A fortunate 
case is the SO (54-43) transition at 206.176062 GHz, which ap- 
pears at 226.572 GHz in the image sideband and permits a con- 
trol of the rejection near the CN N=2-l lines, whose strongest 
hyperfine components are at 226.874 and 226.657 GHz. Single 
sideband system temperatures were in the range 170-190 K in 
December, and 170-250 K in May. 

In December 2011, the exceptionally dry weather and clear 
sky produced rapid cooling a few hours after sunset, which could 
not be compensated by the thermal regulation of the telescope, 
and resulted in significant pointing changes. This also presum- 
ably affected the telescope efficiency. In May, anomalous refrac- 
tion clearly occurred occasionally. As a result, the absolute cali- 
bration is in general accurate to 10-15 % but on several sources 
is not better more than about 25 %. However, since we control 
the receiver sideband rejection, the relative intensities of the de- 
tected lines are expected to be accurate to within 5 % or bet- 
ter, noise being in general a more severe limitation. Given the 
absolute calibration issues, we did not correct for the elevation 
dependent antenna gain, and simply use a uniform conversion 
factor from antenna temperature (TV) to flux density of 9 Jy/K. 

Since our targets should be unresolved by the IRAM 30-m, 
the observations were performed in symmetric wobbler switch- 
ing mode, with the two references +60" away at the same ele- 
vation. Together with the stable weather in December, this mode 
provided flat baselines, allowing clear detection of the contin- 
uum emission from the dust in most sources. Occasional baseline 
discontinuities between the 3 sub-units ("platforming") were 
visible in the spectra, but none of the spectral lines we searched 
for were affected by this problem. 

We observed the sample of 40 sources (see Table [TJ, for 
about 1 hour each. For a typical linewidth of 2.6 km.s~' (2 MHz), 
representative of our sources, this gives a detection sensitivity of 
0.08 Jy. km.s~ (lcr) on the integrated line flux. 

2.3. Data Analysis 

We removed linear baselines in a window ± 60 km.s~' wide 
around each spectral line. Spectra for all objects appear in 
AppendixlBl FiguresQ~]and[2]show spectra for two objects: Haro 
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Table 1. Source sample 



Name 


RA 






Dec 


Spectral 


L, 


5 y ( 1.3 mm) 


Incli. 


A v 


References 




J2000 


.0 




J2000.0 


Type 


Lo 


(mJy) 




O 






FN Tau 


04:14 


14 


59 


28:27:58.1 


M5 


0.50 


< 18 






1.35 


a,A 


CWTau 


04:14 


17 





28:10:56.51 


K3 


0.68 


96 






2.19 


b,B 


CIDA-1 


04:14 


17 


6 


28:06:11 


M5.5 




14 






- 


a,a 


CY Tau (*) 


04:17 


33 


729 


28:20:46.86 


M1.5 


0.40 


111 


28 


±5 


0.1 


c,l,B 


BP Tau 


04:19 


15 


834 


29:06:26.98 


K7 


0.65 


58 


33 


± 6 


0.49 


c,l,B 


DE Tau 


04:21 


55 


6 


27:55:05.55 


Ml 


1.14 


36 






1.56 


b,CJ3 


RY Tau 


04:21 


57 


42 


28:26:35.6 


Kl 


6.59 


229 


66 


±3 


1.84 


b,2,B 


T Tau 


04:21 


59 


435 


19:32:06.36 


K0 


8.91 


200 




30 


1.39 


c,3,C 


Haro 6-5B 


04:22 


01 


00 


26:57:35.5 


K5 


0.05 


134 




80 


9.96 


d,4,DjS 


FT Tau 


04:23 


39 


188 


24:56:14.28 


C 


0.38 


73 


23 


± 5 




c,l,C 


DG Tau-B 


04:27 


02 


56 


26:05:30.4 


c 


>0.02 


531 


64 


± 2 




c,l 


DG Tau 


04:27 


04 


67 


26:06:16.9 


K6 


3.62 


390 


38 


± 2 


3.32 


c,l,D,/? 


Haro 6-10 


04:29 


23 


729 


24:33:01.52 


K5 


2.55 


44 








c,E 


IQ Tau 


04:29 


51 


56 


26:06:44.9 


M0.5 


0.53 


60 






1.25 


a,B 


LkHa 358 


04:31 


36 


15 


18:13:43.1 


K8 


0.09 


17 


28 


±9 


0.13 


a,5,A,y 


HH30 


04:31 


37 


468 


18:12:24.21 


M0 


0.2 - 0.9 


20 


83 


± 2 


2.96 


c,6,F,/? 


HL Tau 


04:31 


38 


413 


18:13:57.55 


K7 


1.53 


819 


45 


± 1 


7.43 


c,l,D,/J 


HK Tau 


04:31 


50 


58 


24:24:17.9 


M0.5 


1.00 


41 




90 


3.41 


b,7,B 


Haro 6-13 


04:32 


15 


419 


24:28:59.47 


M0 


2.11 


114 


40 


± 3 


11.9 


C,5,A,/J 


GG Tau 


04:32 


30 


34 


17:31:40.5 


K7 


0.64 


593 


37 


± 1 


0.76 


b,8,B 


UZ Tau E 


04:32 


43 


071 


25:52:31.07 


Ml 


0.9-1.6 


150 


56 


± 2 


1.49 


c,l,C+G 


04302+2247 


04:33 


16 


2 


22:53:20.0 


- 




130 


88 


± 2 




e,9 


DL Tau 


04:33 


39 


077 


25:20:38.10 


K7 


1.16 


204 


43 


±3 


2.0 


c,l,H 


DM Tau ($) 


04:33 


48 


70 


18:10:10.6 


Ml 


0.16 


109 


35 


± 1 


0.31 


c,l,B,a 


CI Tau (*) 


04:33 


52 


014 


22:50:30.06 


K7 


0.96 


125 


44 


± 3 


1.77 


c,l,B 


AA Tau 


04:34 


55 


42 


24:28:53.1 


K7 


0.66 


73 


70 


± 5 


0.49 


b,10,C 


DN Tau 


04:35 


27 


38 


24:14:58.9 


M0 


0.68 


89 






1.89 


a,D,/? 


DO Tau 


04:38 


28 


59 


26:10:49.5 


M0 


1.29 


136 






2.64 


b,B 


HV Tau C 


04:38 


35 


31 


26:10:38.5 


K6 


0.60 


40 




90 


2.42 


b,ll,B 


LkCa 15 (*) 


04:39 


17 


76 


22:21:03.7 


K5 


0.85 


110 


52 


± 1 


0.62 


c,l,B 


Haro 6-33 


04:41 


38 


827 


25:56:26.68 


M0 


0.76 


34 


52 


± 5 


10.2 


c,5,D 


GO Tau (*) 


04:43 


03 


050 


25:20:18.80 


M0 


0.22 


53 


52 


± 1 


1.18 


a,5,B 


DQ Tau 


04:46 


53 


064 


17:00:00.09 


M0 


0.91 


83 






0.97 


c,C 


DR Tau 


04:47 


06 


22 


16:58:42.8 


K5 


1.97 


159 






2.0 


b,B 


DS Tau 


04:47 


48 


6 


29:25:10.96 


K5 


0.67 


25 






0.31 


b,C 


UYAur 


04:51 


47 


37 


30:47:13.9 


MO 


1.41 


29 




42 


2.05 


b,12,B 


GM Aur (*) 


04:55 


10 


98 


30:21:59.5 


K7 


1.23 


176 


50 


± 1 


0.14 


c,l,B 


AB Aur 


04:55 


45 


80 


30:33:04.0 


B9 


44.86 


110 


2C 


-30 


0.62 


f,13,I 


SUAur 


04:55 


59 


4 


30:34:01.39 


G2 


9.29 


< 30 




62 


0.90 


b,14,B 


MWC 480 


04:58 


46 


27 


29:50:37.0 


A2 


11.50 


289 


37 


± 1 


0.25 


c,l,G,y 


VDA98 CB 26 


04:59 


50 


74 


52:04:43.8 






190 




90 




g,15 


CIDA-8 


05:04 


41 


4 


25:09:57 


M3.5 


0.26 


8 






0.83 


a,A,a 


CIDA-11 


05:06 


23 


3 


24:32:24 


M3.5 


0.22 


< 4 






0.36 


a,A 


RWAur 


05:07 


49 


56 


30:24:05.1 


K3 


1.72 


42 


45 


±5 


0.32 


b,16,B 


MWC 758 


05:30 


27 


51 


25:19:58.4 


A3 


11.00 


55 




21 




c,17,J 


CQ Tau 


05:35 


58 


485 


24:44:54.19 


A8/F2 


12.00 


162 


29 


± 2 


0.96 


c,l,J,y 


HD 163296 (#) 


17:56 


21 


31 


-21:57:22.0 


Al 


37.70 


670 


45 


±4 


0.25 


h,K,y 



Notes. Data for sou rces marked with (*) are for CN only, and were obtain ed s eparately, ($) no SO lin e, ( #) different SO line. 

Flux referenc es: (a)ISchaefer et al.l ((2009) (b) A ndrews & Williams! J2001) (c)lGuilloteau et alj d2011l) (d) IOsterloh & Beckwithl d 19951) (e) Graeffe 

et al 2012, (ft lPietu et al] d2005h. (g)lLaunhardt et al] d2009h (h)lQi et alJ 1201 

Inc lination reference s: (l )lGuilloteau et alJfeOl ll T( 3)llsella et alld2010bl) ( 3) lRatzka etai1d2009l) . (4) lKrist et alldl998l). (5)ISchaefer et all (2 009) 
(6)lPetv et all d2006h. (7)|McCabe et alj d201 lh. (8) Guillot eau et al .Id 19991). (9) Gra effe et al 2012, (lO)lBouvier et alj dl999l). (1 llDuchene et al] 
<201Ch. (12) lDuvert et all d!998l). (13)lPietu et al] d2005l) . (14) lAkeson et al] ( [20021) (15) lLaunhardt et al] d2009l) . (i6) ICabrit et al] d2006l) . (17) 
llsella et al] d2010an . (18) llsella et al] d2007h . 

Lu minosity reference s: (A ) ISchaefer et alj d2009T). (B) iBertout et alj d2007l). (C)lKenvon & Hartmannl d 19951). (D)IWhite & Hillenbrandl d2004l) 
(E)lPrato et all d2009h. (F)lPety et al] d2006h . (G) ISimon et al.l d2000h . (H) IWhite & Ghezl d200lh . (I) lHernandez et alj d2004h . (J) IChapillon et al] 
d2008l) . (K) lTilling et al] d2012l) 

Extinction references : iKenvon & Hartmannl dl995h . except (a) iBriceno et al] d2002h (Aj for CIDA-8), (fi) IWhite & Hillenbrandl d2004h . (y) 
Ivan den Ancker et all (1998) 



6-10 which exhibits bright lines in all molecules, and MWC 480, 
an isolated HAe s tar for which CN has a lready been imaged with 
the IRAM array (Chapi llon et al.ll2~012l) . included in our sample 
as a reference object. To derive line intensities, we fitted simple 



Gaussian profiles to the observed lines taking into account the 
hyperfine structure of the C I7 J=2-l and CN N=2-l lines. We 
used the "HFS (HyperFine Structure)" method of the CLASS 
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-20 20 

Velocity (km/s) 

Fig. 1. Spectra of the observed transitions towards Haro 6-10. 
Top panel: 13 CO J=2-l spectrum. Bottom panels, from top to 
bottom: spectra of 13 CO J=2-l, C n O J=2-l, CN N=2-l, ortho 
H 2 CO 3 13 -2i2, ortho H 2 CO 3 l2 - 2 n , and SO 5 4 -4 3 presented 
on a common scale. The Y axis scale is antenna temperature, T^ 
(the conversion factor to flux density is 9 Jy/K). 




20 
Velocity (km/s) 



Fig. 2. As FigQ]but towards MWC 480. The red curve is a simple 
Gaussian fit to the disk emission in I3 CO J=2-l. 



data analysis program in the GILDASQ package for this purpose. 
This method assumes a common excitation temperature for all 
the hyperfine components, a likely situation at the high densities 
appearing in disks. For a Gaussian line, the emerging profile is 
then 



7» = /7(^CT ex )-./ v (7bg)): 



(1) 



2/ 



1 - exp 



—RfT exp 



/ V - V sys - V; 

SV 



2\\ 



where are the relative intensities = 1), V; the velocity of 
each hyperfine component with respect to a common reference 
frequency, and t the total line opacity (sum over all hyperfine 
components, as the R, are normalized). V sys is the source veloc- 
ity, 5V the line width at 1/e, and r\ a beam filling factor. The 
HFS method uses X = Tjr(J Y (T ex ) - J v (T\, g ), V sys , 6V and r as 
free parameters. The difference in saturation between the vari- 
ous hyperfine components allows determination of r provided 
sufficient signal to noise. For CN, although all 19 hyperfine 
transitions were used for the fitting, the quoted integrated line 
flux includes only the 3 blended hyperfine components around 
226.874 GHz and the two weak components at 226887.399 & 
226892.151 MHz (which appear near V LSR -11 and -17 km.s 1 
in the brighter sources) Q. Given the line ratios, the total line flux 
including all hyperfine components would be 1.67 times larger 
in the optically thin limit. 



Table 2. Frequencies of transitions 



Frequency 
(MHz) 



Molecule and 
Transition 



Detection 



206176.062 








SO 


(5 4 


-4 3 


) 




Yes 


209230.201 






HC 3 N J= 


23-22 




No 


209419.15 








I3 C 


8 o 


r=2- 


1 




cloud only 


211211.455 




ortho H 


,co 


3l3 




In 


Yes 


220398.688 








13 CO J 


=2-1 






Yes 


224714.385 








c r 


OJ 


=2-1 






Yes 


225697.781 




ortho H 


,co 


3n 




In 


Yes 


226287.4185 


CN 


N= 


=2- 


1 J= 


3/2 


3/2, 


F= 


=1/2-1/2 




226298.9427 


CN 


N= 


=2- 


1 J= 


3/2 


3/2, 


F= 


=1/2-3/2 




226303.0372 


CN 


N= 


=2- 


1 J= 


3/2 


3/2, 


F= 


=3/2-1/2 




226314.5400 


CN 


N = 


=2- 


1 J= 


3/2 


3/2, 


F= 


=3/2-3/2 


Yes 


226332.4986 


CN 


N= 


=2- 


1 J= 


3/2 


3/2, 


F= 


=3/2-5/2 




226341.9298 


CN 


N = 


=2- 


1 J= 


3/2 


3/2, 


F= 


=5/2-3/2 




226359.8710 


CN 


N = 


=2- 


1 J= 


3/2 


3/2, 


F= 


=5/2-5/2 


Yes 


226616.5714 


CN 


N = 


=2- 


1 J= 


3/2 


1/2, 


F= 


=1/2-3/2 


Yes 


226632.1901 


CN 


N= 


=2- 


1 J= 


3/2 


1/2, 


F= 


=3/2-3/2 


Yes 


226659.5584 


CN 


N= 


=2- 


1 J= 


3/2 


1/2, 


F= 


=5/2-3/2 


Yes 


226663.6928 


CN 


N = 


=2- 


1 J= 


3/2 


1/2, 


F= 


=1/2-1/2 


Yes 


226679.3114 


CN 


N = 


=2- 


1 J= 


3/2 


1/2, 


F= 


=3/2-1/2 


Yes 


226874.1908 


CN 


N= 


=2- 


1 J= 


5/2 


3/2, 


F= 


=5/2-3/2 


Yes 


226874.7813 


CN 


N= 


=2- 


1 J= 


5/2 


3/2, 


F= 


=7/2-5/2 


Yes 


226875.8960 


CN 


N = 


=2- 


1 J= 


5/2 


3/2, 


F= 


=3/2-1/2 


Yes 


226887.4202 


CN 


N = 


=2- 


1 J= 


5/2 


3/2, 


F= 


=3/2-3/2 


Yes 


226892.1280 


CN 


N= 


=2- 


1 J= 


5/2 


3/2, 


F= 


=5/2-5/2 


Yes 


226905.3574 


CN 


N= 


=2- 


1 J= 


5/2 


3/2, 


F= 


=3/2-5/2 




227191.8195 


CN 


N= 


=2- 


1 J= 


5/2 


1/2, 


F= 


=5/2-3/2 





Note s. CN N=2-l l ine frequencies were measured in laboratory by 
ISkatrud et all dl983h; we us e here the fitted values from the CDMS 
Database dMiiller et alJ2001h . 



1 See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about GILDAS softwares 

2 This is essentially the flux in the J=5/2-3/2 fine structure line, as the 
only missing hyperfine component has much smaller intensity 
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0.1 1 _ 10 

o-H 2 C0 3 12 ™2 11 line flux (Jy.km.s ) 



Fig. 3. Integrated line flux of the 3n - 2n line of H2CO vs that 
of the 3 13 - 2n transition. Color codes indicate detected vs non- 
detected lines: blue for detection of both lines, black for detec- 
tion of the 3 13 - 2 12 line, green for non detections. The detection 
threshold is set at 4<x. For detected lines, the errorbars are lcr. 
For non-detected lines, the limit is the 2<x value, and the arrow 
length is lcr. Sources in italics are those with known molecular 
outflows. 



The two transitions of ortho-F^CO exhibit a strong correla- 
tion (see Figf3]l because they arise in the same rotational levels 
and differ only by the small K splitting, with the 21 1.2 GHz line 
on average 20% stronger than the 225.2 GHz line. We therefore 
computed an average spectrum and line intensity using this ra- 
tio. For the several lines we analyzed, Tables [3]and|4]list the de- 
rived line intensity, Wlsr and line width, AV. Figures |4]through 
Qdisplay relationships of the line intensities. The color coding in 
these plots is as follows: blue designates detection of both lines, 
red designates detection only of the line plotted on the x-axis, 
and green detection only of the line plotted along the y-axis, and 
black is for non detection of both lines. Sources in italics are 
those with known molecular outflows. 

The 13 CO results are seriously affected by confusion. In 
many cases, we attempted to fit only a putative disk component, 
by masking the channels in which confusion seemed likely, ei- 
ther because they appear negative or exceedingly strong, or be- 
cause of strong C n O emission at the same velocity. A summary 
of disk detection is given in Table The fits are displayed on 
top of the 13 CO spectra in the Figures of Appendix IE1 The fit 
results can be found in Table [3] These derived line intensities 
are presented as a function of CN line flux in Fig |9] Because of 
the masking process to ignore confused parts of the spectrum, 
the errorbars can be large. However, in sources where both 13 CO 
and CN are detected, they in general display similar intensities, 
within a factor 2. 

3. Results 

Line emission from these objects can come a priori from 4 
distinct regions: molecular cloud(s) along the line of sight, the 



circumstellar disk, a molecular outflow and (specially for the 
youngest objects) a remnant envelope. Molecular cloud will ex- 
hibit narrow lines, and because of our observing technique which 
differentiate the emission with that 1' away, may appear in emis- 
sion or absorption. Disks should lead to symmetric, in general 
double-peaked, line profiles, and line widths oc sin i (with typi- 
cal values 2-4 km.s -1 , face-on objects being rare). Contribution 
from outflows or envelopes may vary considerably, and may be 
difficult to distinguish from disk emission. Most of the stars in 
Table Q] have a rich observational history from optical, infrared 
and radio molecular line observations which can help interpret- 
ing the line profiles. We have compiled the relevant informa- 
tion and comment on the observational results for each source in 
Appendix lAl 



3.1. Confusion in 13 CO 

The 13 CO spectra are in general complex. Out of the 42 sources, 
the only two in which 13 CO exhibits the expected double-peaked 
profile which is a clear signature of emission from Keplerian 
disk are the previously studied objects GG Tau and MWC 480 
(see FigJ3J, included here as reference objects. CQ Tau and 
MWC 758 also exhibit 13 CO spectra which can be attributed 
solely to disk emission. These 4 sources had been previously 
studied specifically because they are relatively far from the 
molecular clouds. 

For the other sources, the spectrum complexity has two main 
origins. Narrow features appear in emission or apparent absorp- 
tion. These narrow features result from small velocity and in- 
tensity gradients in the molecular cloud which are on the same 
line of sight. Our beam-switching observing procedure reveals 
such gradients on the arcminute scale. Note that, as the switch- 
ing was made at constant elevation, the gradient direction rotates 
on the plane of the sky, so that the observed features actually de- 
pend on when the observations were performed. A comparison 
bet ween the spe c tra of BP Tau obtained here and that published 
by iDutrev et al.l d2003l) illustrates this effect. The BP Tau case 
also demonstrates that the lack of apparent confusion does not 
imply the source is unaffected by the cloud: a sufficiently opti- 
cally thick and homogeneous cloud will hide any disk emission 
situated behind. 

Besides contamination by molecular clouds, many sources 
also exhibit strong, relatively wide profiles in 13 CO. Based on 
prior knowledge about these sources, this emission is in gen- 
eral coming from outflows or envelopes, as in HL Tau, T Tau, 
Haro 6-10, CB26, DG Tau and DG Tau B, which are all known 
to drive powerful jets and CO outflows (see Appendix|A]i. This 
interpretation is corroborated by the emission lines from other 
molecules (see below). However, the line widths are not very 
large, so this emission can also contain a substantial contribu- 
tion from the disk, specially for edge-on disks. An example is 
IRAS04302+224 7 (hereafter 04302+2247, the Butterfly star, see 
IWolf et aLll2003h . for which IRAM interferometric observations 
revea l that the detected lines clearly trace the disk (Dutre vet al.l 
I2013L in prep). 

Despite confusion with molecular clouds, many other 
sources exhibit line wings which are potential signatures of disk 
emission. As mentioned in Sec l2.31 these sources are identified 
in Table|5] and appear in Fig|9] 
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Fig. 4. Integrated line flux of H2CO (combina- 
tion of both transitions) vs CN N=2-l (main 
group of hyperfine components). Color codes 
indicate detected vs non-detected lines: blue for 
detection of both lines, red for detection of the 
x-axis line only (here CN), black for detection 
of the y-axis line only (here H 2 CO), and green 
for non detections. The detection threshold is 
set at 3cr. For detected lines, the errorbars are 
lcr. For non-detected lines, the limit is the 2a 
value (best fit + 2<r), and the arrow length is lcr. 
Sources in italics are those with known molec- 
ular outflows. 



3.2. CN as a confusion-free disk tracer 

There is essentially no evidence for confusion with molecular 
clouds in CN, except for a few (4 or 5) sources. This can be seen 
by the lack of narrow features corresponding to those detected 
in 13 CO. 

Only two sources exhibit clear evidence of confusion in CN: 
CW Tau and CIDA-1 (see Fig|R2] and |R3]l in which the CN 
lines appear in absorption. These two sources are separated by 
a few arcmin only, and all other lines also display the same be- 
haviour, clearly demonstrating the confusion along this line of 
sight. The lines are very narrow, 0.44 + 0.04 km.s~', a further 
indication that this comes from a cold molecular cloud. More 
detailed arguments about the properties of the molecular cloud 
and its impact as a confusing signal for disk emission are given 
in AppendixICl 

A closer examination reveals confusion in two other sources, 
AA Tau and AB Aur, and an unclear situation for Haro 6-13, 
which exhibit emission with relatively narrow lines (see Table 
[3] and Appendix lAl for details). Finally, like for 13 CO, a lack of 
apparent confusion does not indicate an absence of confusion. 

CN emission is compatible with optically thin emission in 
most sources (r p < 0.5, where r p is the peak line opacity). 
The HFS method yields r, the sum of opacities of all hyperfine 
components. The main group of components accounts for a lit- 
tle more than half of this sum, but with local linewidths in the 
range 0.2-0.5 km.s -1 , the peak opacity is t p ~ 0.27 - 0.4t, be- 
cause the hyperfine components overlap only very moderately, 
being separated by 0.77 and 1.50 km.s -1 . The opacity is con- 
strained in the brightest sources: T Tau has optically thin emis- 
sion (r = 0.53 ± 0.06) but Haro 6-10 exhibits substantial optical 
depth (t = 3.0 + 0.3, i.e. t p ~ 1). The next brightest source, DG 
Tau, is marginally thick (r = 0.9 + 0.3). The spectra also suggest 



CN could be optically thick in the Butterfly star 04302+2247, 
but because of limited signal-to-noise ratio, the significance is 
low (r = 10+3). Spectra of these sources are displayed in FigfTUl 
together with the best fit profiles. 

Non zero optical depth (r ^ 1) is also suggested for some of 
the other strong sources, such has DL Tau, and also DM Tau 
and LkCa 15, but the err or bars remain l arge. The interfero- 
metric study performed by Chapil lon et alj d2012l) however con- 
firms that CN is not completely optically thin in the latter two 
sources. Apart from Haro 6-10 and 04302+2247, we found no 
other source in our sample with optical depth larger than ~ 1 ; 
the weighted mean opacity is ~ 0.60- 

In summary, we detected CN N=2-l in 24 out of the 42 
stars in our sample, and have only 4 cases of confusion, and 
essentially all spectra are consistent with total optical depth 
-0.3-0.6. 

3.3. SO and H 2 CO as outflow or envelope tracers 

Strong SO emission is found in 7 out of the 41 stars in which 
it has been observed (see Fig|6]l. All these sources also display 
relatively strong H2CO emission (see FigQ). All of them, except 
04302+2247, are driving well-known jets and outflows. In these 
sources, CN is either quite strong (DG Tau, T Tau, DG Tau B, 
Haro 6-10) or very weak (CB26, HL Tau). All these sources are 
also deeply embedded, and surrounded by an extended complex 
envelope. We also have a marginal detection of SO in Haro 6- 
5 B, another embedded, jet-driving object. The line profiles of 
SO and H2CO are similar, and resemble those of C 17 when 
detected. As the line widths are not very wide, this emission may 



3 The same value is obtained with or without T Tau. 
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CN N = 2-1 J = 5/2-3/2 line flux (Jy.km.s ) 



Fig. 5. As for FigHbut for 13 CO J=2-l (includ- 
ing contribution from the cloud) vs CN N=2-l. 
T Tau is out of scale here. 




0.1 1 10_ 

CN N = 2-1 J=5/2-3/2 line flux (Jy.km.s 1 ) 



Fig. 6. As for FigHbut for SO vs CN N=2-l. 
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be attributed to the envelope rather than the outflow itself. On the 
contrary, CN lines often display profiles which are quite different 
from those of SO, being either narrower (e.g. HL Tau, Fig JB. 16t 
or with no narrow velocity component (DG Tau B, Fig lB.lOb . or 
(rarely) broader (DG Tau, see Fig JB. lib . 

The detection rate of H2CO is rather high, with 20 out of 
42 stars observed (not counting the dubious cases FN Tau and 
CIDA-8). However, cloud contamination appears more impor- 
tant in H2CO (7 sources with narrow lines out of 42) than in 
CN. This can be explained by excitation considerations, because 
of the smaller critical densities of the observed transitions of 
H2CO compared to CN N=2-l. These transitions have similar 
Einst ein coefficients, but th e collisional rates of H2CO with para- 
H2 dTroscompt et al.l2009h are typ ically 5 to 10 times l arger than 
those of CN with para-H 2 at 10 K dKalugina et al.ll2012l) . On av- 
erage, CN is brighter than H2CO (see FiglUi, and the high detec- 
tion rate of H2CO is obtained because two lines of similar inten- 
sities could be observed simultaneously and added in this study. 
Strong H2CO lines clearly correlates with SO detectability (see 
FigQ), i.e. in general outflow-driving sources. 

3.4. C 17 O and other lines 

While 13 CO emission from disk is overwhelmed by confusion, 
C 17 is in general too faint to be a suitable tracer of disks 
(only 9 detections out of 40 observed sources, see Fig|8]l. It is 
convincingly detected in the disks of AB Aur, MWC 480 and 
04302+2227 only. On another hand, residual confusion from 
the cloud is still visible in many sources: CW Tau and CIDA- 
1, CIDA-8, Haro 6-13, HK Tau, Haro 6-33, Haro 6-5 B, HH 30, 
and LkHa 358. The strong emission from HL Tau and CB26 pre- 
sumably comes from their outflows, as the lineshapes are com- 



plex and resemble those of SO. For other "outflow" sources, the 
linewidths are small enough to leave some doubt about the ori- 
gin. 

We also have a 5<x detection of 13 C 18 towards Haro 6-10, 
with a relatively narrow linewidth of 0.7 ±0.15 km.s . A line of 
similar width is also detected towards AB Aur, at the Act level. In 
the later source, as the disk emission spreads over 2-3 km.s 1 , 
this must come from the molecular cloud, and implies that this 
cloud also contributes to the signal of the other isotopologues. 

Finally, emission at the 3<x level appears towards Haro 6-10, 
MWC 480 and GG Tau near the expected frequency of the HC3N 
J=23-22 transition, but these features are either too broad, too 
narrow, or not at the right velocity. No other significant feature 
appears anywhere else in the observed bandwidth. 

4. Discussion 

4. 1 . Using CN as a disk tracer 

Our main approach to test whether CN is a disk tracer is to as- 
sume that the (mean) CN surface density is independent of the 
disk properties. This assumption is sustained by several inde- 
pendent arguments. Firs t, the imaging s tudy o f DM Tau, LkCa 
15 and MWC 480 by IChapillon et atl (1201 2h give similar CN 
surface densities in the 3 objects. Second, chemical modeling 
performed to interpret these CN images predict that the CN sur- 
face density is rather insensitive to the disk mass, varying by 
only ±30% for masses differi ng a factor 3 i n bot h directions. A 
similar result was found by Jonkheid et al. (2007) in their mod- 
els for disk around HAe stars. Last, from an observational point 
of view, our analysis shows that the CN line optical depth does 
not appear to exhibit wide variations. We thus assume the CN 
line opacity to be constant (and relatively small), apart from the 
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Table 3. Fit Results: CN and 13 CO 



Source CN "CO 

Area V LSR AV Area V LSR AV 

Jy. km.s~' km.s~' km.s 1 Jy.km.s~ 1 km.s 1 km.s~' 



FNTau 




< 





38 


















-1.40 


¥0.09 


7.44+ 0.03 


0.84 


+ 0.05 


Cloud 


CWTau 


-0.66 


± 





07 


6 


80 = 


= 


07 





85 = 


= 


07 


2.69 


+ 0.45 


6.20 


±0.10 


3.25 


±0.38 




CIDA-1 


-0.82 


± 





09 


6 


93 = 


= 


02 





42=i 


= 


03 


1.14 


±0.75 


8.23 


±0.19 


0.56 


±0.27 




BPTau 




< 





59 


















0.26 : 


t 0.14 


3.65± 0.39 


1.06 


+ 0.60 


Cloud 


DETau 




< 





39 


















0.20 


+ 0.15 


11.43 


±0.99 




[ 3.00] 




RYTau 


3.91 


± 





86 


5 


06 = 


= 1 


58 


14 


22=i 


= 2 


19 


2.54 


± 0.51 


7.24 


±0.29 




[ 3.00] 




TTau 


18.62 


± 





31 


7 


96 = 


= 


03 


3 


12=i 


= 


05 


170.37-. 


t 0.20 


8.02± 0.00 


2.50 


±0.00 


Cloud 


Haro 6-5B 


0.61 


± 





13 


7 


67 = 


= 


42 


3 


61 =i 


= 


73 


14.94 : 


± 0.30 


7.88± 0.01 


0.80 


±0.02 


Cloud 


FTTau 


0.63 


± 





12 


6 


95 = 


= 


28 


2 


21 =i 


= 


42 


0.15 


+ 0.07 


9.38 


+ 0.20 


0.78 


±0.34 




DG Tau-B 


1.33 


± 





15 


6 


18 = 


= 


18 


2 


34=i 


= 


22 


40.59 - 


t 0.11 


6.43± 0.00 


2.03 


±0.00 


Cloud 


DGTau 


6.62 


± 





16 


6 


52 = 


= 


03 


2 


28=i 


= 


07 


42.08-. 


± 0.14 


6.67: 


t 0.00 


2.20 


±0.01 


Cloud 


Haro 6 


11.05 


± 





22 


6 


70 = 


= 


02 


1 


97 =i 


= 


05 


41.64 : 


k0.47 


7.15± 0.01 


1.07 


±0.02 


Cloud 


IQTau 


2.28 


± 





21 


6 


22 = 


= 


25 


5 


07 =i 


= 


35 


0.45 


+ 0.14 


4.70 


+ 0.24 


1.45 


±0.60 




LkHa 358 




< 





27 


















6.91 


± 1.31 


6.82 


±0.04 


1.43 


±0.12 




HH 30 




< 





41 


















8.21 


± 2.97 


7.18 


±0.07 


2.34 


±0.29 




HLTau 


0.81 


± 





13 


5 


00 = 


= 


20 


1 


91 =i 


= 


37 


54.88-. 


t 0.76 


5.88± 0.01 


1.46 


±0.03 


Cloud 


HKTau 


1.33 


± 





22 


4 


49 = 


= 


76 


9 


25 =i 


= 1 


14 


1.40 


+ 0.62 


5.25 


+ 0.62 


3.12 


±0.77 




Haro 6-13 


0.30 


± 





08 


5 


37 = 


= 


14 





69=i 


= 


16 


10.05 


± 7.19 


5.03 


±0.06 


1.67 


±0.12 




GGTau 


3.20 


± 





34 


6 


31 = 


= 


17 


2 


40d 


= 


22 


5.01 


±0.34 


6.49 


±0.09 


2.66 


±0.17 




UZ Tau E 


1.21 


± 





19 


5 


96 = 


= 


47 


5 


75 =i 


= 


66 


2.61 


±0.29 


5.90 


±0.18 


4.47 


±0.46 




04302+2247 


1.53 


± 





24 


5 


92 = 


= 


48 


5 


37 =i 


= 


55 


33.71 


+ 1.78 


5.62 


±0.02 


2.92 


±0.06 




DL Tau 


1.82 


± 





15 


6 


70 = 


= 


10 


1 


41 =i 


= 


18 


0.53 


± 0.24 


6.62 


±0.29 


1.29 


±0.64 




DM Tau 


5.90 


± 





34 


6 


04 = 


= 


07 


1 


86=i 


= 


12 


6.84 


± 0.17 


6.07 


±0.02 


1.58 


±0.04 




AA Tau 


0.74 


± 





10 


7 


36 = 


= 


11 





96 i 


= 


15 


1.72 


±0.28 


6.38 


±0.53 


6.16 


± 1.03 




DN Tau 


1.54 


± 





17 


6 


61 = 


= 


20 


3 


17 =i 


= 


29 


0.24 


+ 0.09 


7.79 


±0.22 


0.98 


±0.44 




DO Tau 


0.84 


± 





20 


5 


91 = 


= 


22 


1 


44d 


= 


54 


1.18 


± 0.35 


3.61 


±0.85 


4.96 


± 1.65 




HV Tau C 


1.22 


± 





20 


6 


95 = 


= 


59 


6 


51 =i 


= 


67 


2.14 


± 0.61 


6.08 


±0.34 


5.02 


±0.87 




Haro 6-33 


0.91 


± 





17 


5 


88 = 


= 


33 


2 


98 i 


= 


53 


1.59 


±0.60 


4.95 


±0.55 


4.52 


±0.98 




DQ Tau 




< 





34 


















-13.78-. 


t 0.14 


9.62± 0.01 


1.62 


±0.01 


Cloud 


DR Tau 




< 





75 


















6.73 ± 


10.82 


9.79 


± 1.95 


1.80 


± 1.42 




DS Tau 




< 





28 


















0.13 -. 


t 0.08 


8.13± 0.22 


0.64 


±0.45 


Cloud 


UY Aur 




< 





69 


















24.09 


+ 0.48 


6.30 


±0.01 


1.41 


± 0.02 




AB Aur 


0.42 


± 





13 


6 


87 = 


= 


08 





47 =i 


= 


09 


36.43 


+ 0.17 


5.89 


±0.00 


1.56 


±0.01 




SUAur 




< 





27 


















1.75 


±0.36 


7.09 


±0.33 


3.05 


±0.47 




MWC 480 


3.58 


± 





23 


5 


04 = 


= 


16 


4 


44=i 


= 


22 


7.42 


± 0.21 


5.04 


± 0.06 


4.18 


± 0.13 




CB26 


0.82 


± 





14 


6 


20 = 


= 


35 


3 


35=i 


= 


45 


48.36 -. 


t 0.15 


5.6<S± 0.00 


1.39 


±0.00 


Cloud 


CIDA-8 




< 





26 


















0.46 


+ 0.23 


7.69 


±0.68 




[ 3.00] 




CIDA-11 




< 





23 


















0.11 -. 


t 0.05 


9.04± 0.11 


0.30 


±0.33 


Cloud 


RWAui 




< 





60 


















0.84 


± 0.31 


2.46 


±0.59 


2.86 


± 1.70 




MWC 758 




< 





86 


















0.82 


±0.25 


6.27 


±0.39 


2.28 


±0.61 




CQ Tau 




< 





28 


















1.54 


±0.18 


6.01 


±0.34 


5.65 


±0.70 




HD 163296 


10.91 


± 





23 


5 


85 = 


= 


05 


4 


28 =i 


= 


07 


6.87 


±0.68 


6.05 


+ 0.20 


4.02 


±0.45 





Notes. Integrated line flux of CN (main group of hyperfine components) and 13 CO. For 13 CO, results in italics include cloud (or outflow) emission, 
while the others were obtained by masking confused channels and trace the disk emission only. 



inversely proportional to cos(0, so the dependency on inclina- 
tion disappears, and the flux just scales with the total number 
of molecules. In this optically thin case, Eq|2] can be extrapo- 
lated to low inclinations i < 20°. For nearly edge-on objects , 
where the above formula breaks, iBeckwith & Sargenl <[l993) 
have shown that the expected flux is not a strong function of 
inclination. Hence the apparent dependence on i essentially van- 
ishes. We thus use pcos(0 = 0.3cos(45°) for all sources in 
which the opacity cannot be determined from our spectra, and, 
for the few sources in which r is determined by the HFS fit, 
we use pcos(/) = rcos(45°). However, larger disks with smaller 
average line opacity can equally fit the data: in this respect, the 
derived radii are thus likely to be lower limits to the true CN disk 
size. 

For the other parameters, we take AV = 0.2km.s _1 , a value 
consistent with the local line width determined in the disks 



3 sources in which it is clearly determined. Under this assump- 
tion, the total flux is directly related to the total number of CN 
molecules, i.e. principally by the disk outer radius (S v cc Rg Ut ). 

For a quantitative derivatio n of R out , we use the approach 
of Guilloteau & Dutrevl d!998l) which shows that, at moderate 
inclinations (20° < i < 70°), the line flux from a Keplerian disk 
is related to local line width and disk outer radius: 



Rcmt — D 



J Sydv 



fl v (7o)(pAV)jrcos(i) J 



1/2 



(2) 



where D is the distance (140 pc), Tq the excitation tempera- 
ture, AV the local line width, i the inclination, and p a factor 
depending on line opacity along the line of sight. The rationale 
for and derivation of Eq|2]are given in Appendix IPJ1 In the opti- 
cally thin regime, p is proportional to the opacity, which is itself 
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Table 4. Fit Results: H 2 CO and C I7 



Source 




H 2 CO 








C 


n „ ■ 

"O 








Area 




Vlsr 




AV 


Area 




Vlsr 




AV 




Jy. km.s~' 




fcm.s -1 




fcm.s -1 


Jy. km.s~' 




fcm.s -1 




km.s~' 


FNTau 


0.24 ± 0.06 


-0.53 


± 0.14 


1.10 


± 0.31 


< 0.20 










CWTau 


-0.27 ± 0.03 


6.80 


±0.03 


0.51 


±0.05 


0.21 ±0.05 


6.28 


±0.06 


0.43 


±0.11 


CIDA-1 


-0.25 ± 0.04 


6.98 


±0.03 


0.43 


±0.09 


-0.87 ± 0.07 


6.91 


±0.02 


0.43 


±0.01 


BPTau 


<0.36 










< 0.34 










DETau 


<0.26 










<0.56 










RYTau 


< 0.28 










< 1.68 










TTau 


12.68 ± 0.14 


7.87 


±0.01 


2.51 


± 0.03 


6.07 ± 0.27 


7.96 


±0.04 


2.03 


±0.11 


Haro 6-5B 


0.19 ±0.04 


8.40 


± 0.21 


1.49 


± 0.29 


< 0.84 










FTTau 


< 0.25 










< 1.08 










DG Tau-B 


1.54 ±0.12 


6.43 


± 0.05 


1.61 


± 0.19 


4.01 ± 0.21 


6.44 


±0.02 


0.98 


±0.06 


DGTau 


2.37 ± 0.08 


6.43 


±0.03 


1.50 


±0.06 


1.34 ±0.22 


5.87 


±0.15 


2.24 


±0.48 


Haro 6 


7.70 ± 0.09 


6.59 


±0.01 


1.44 


±0.02 


4.39 ±0.18 


6.39 


±0.02 


0.98 


±0.05 


IQTau 


< 0.25 










< 0.23 










LkHa 358 


0.12 ±0.03 


7.14 


± 0.03 


0.28 


± 0.57 


1.09 ± 0.10 


7.05 


± 0.02 


0.45 


±0.04 


HH 30 


0.18 ±0.04 


6.06 


± 0.11 


0.78 


± 0.15 


1.81 ±0.06 


6.52 


±0.01 


0.43 


±0.00 


HLTau 


2.25 ± 0.12 


6.66 


± 0.12 


4.25 


± 0.28 


5.27 ± 0.23 


6.35 


± 0.11 


5.08 


±0.28 


HKTau 


< 0.21 










< 0.18 










Haro 6-13 


1.13 ± 0.10 


5.81 


± 0.17 


3.85 


±0.43 


2.63 ± 0.30 


5.60 


±0.01 


0.43 


±0.01 


GG Tau 


0.91 ± 0.19 


6.40 


± 0.24 


2.02 


±0.46 


< 0.50 










UZ Tau E 


<0.25 










< 1.11 










04302+2247 


2.82 ±0.15 


5.73 


±0.13 


4.88 


± 0.27 


2.23 ± 0.24 


6.24 


± 0.20 


3.77 


±0.51 


DL Tau 


0.65 ±0.11 


6.59 


±0.26 


2.98 


±0.54 


< 0.21 










DM Tau 


0.54 ± 0.05 


6.21 


±0.06 


1.22 


± 0.13 












AA Tau 


< 0.46 










< 0.22 










DN Tau 


0.26 ± 0.08 


5.73 


±0.28 


1.62 


± 0.56 


< 0.21 










DO Tau 


0.33 ±0.11 


5.35 


±0.39 


2.15 


±0.54 


<0.68 










HV Tau C 


< 0.18 










< 0.31 










Haro 6-33 


0.36 ± 0.05 


5.82 


±0.05 


0.67 


± 0.12 


< 0.24 










DQ Tau 


< 0.20 










<0.46 










DR Tau 


< 1.36 










< 1.61 










DS Tau 


0.49 ±0.12 


6.08 


± 0.44 


3.46 


±0.80 


< 0.22 










UY Aur 


0.78 ± 0.17 


5.92 


± 0.21 


1.79 


±0.40 


<0.50 










AB Aur 


0.90 ±0.12 


6.11 


±0.19 


2.56 


± 0.29 


2.61 ± 0.28 


6.02 


±0.10 


2.03 


±0.25 


SU Aur 


< 0.16 










< 0.26 










MWC 480 


< 0.14 










1.66 ± 0.26 


5.51 


±0.46 


5.60 


±0.88 


CB26 


1.56 ±0.09 


5.78 


±0.06 


2.45 


± 0.20 


5.77 ± 1.39 


5.59 


±0.04 


0.82 


±0.11 


CIDA-8 


0.16 ±0.03 


6.75 


±0.02 


0.28 


± 0.58 


0.77 ± 0.07 


6.77 


±0.02 


0.43 


±0.00 


CIDA-11 


< 0.19 










< 0.18 










RW Aur 


< 0.43 










< 0.73 










MWC 758 


< 0.51 










2.51 ± 0.77 


-0.31 


± 3.07 


17.47 


±4.84 


CQ Tau 


0.51 ± 0.10 


6.15 


±0.49 


4.41 


± 0.77 


< 2.47 










HD163296 


1.38 ±0.20 


5.50 


± 0.37 


5.12 


± 0.79 













Notes. Integrated line flux of H 2 CO (average of both transitions) and C 17 



of DM Tau, LkCa 15 an d MWC 480 by iPietu et alJ (120071) 
and Chapillon et alJ (12.0121) . For the temperature, Chapillon et al.l 
(120121) derived 30 K for MWC 480, but values (10 - 15 K) for 
DM Tau and LkCa 15. We used 15 or 30 K according to our best 
knowledge of source properties. Note that the derived disk ra- 
dius will scale as 1/ V^o (Eq|2]). The assumed temperatures and 
derived disk radii are given in Table[5]for a uniform disk distance 
of 140 pc. Given the uncertainty on the CN opacity, these radii 
should only be interpreted as estimates of the disk size. They 
range from < 200 AU for non detection, to > 1000 AU for the 
brightest sources. For MWC 480 and DM Tau, the derived ra- 
dius is in good agreement wi th the more accurate d etermination 
from interferometric data of Cha pillon et alJ (1201 2l) . justifying a 
posteriori our choice of parameter values. 

For T Tau and Haro 6-10, the CN line opacity is reason- 
ably well constrained. The large disk radii which are derived 
make the disk interpretation very unlikely, as both objects con- 



tain multiple systems with strong evidence for tidal truncation 
of the embedded circumstellar disks to ou ter radii < 60 AU 
dAkeson et alJll998l: iGuilloteau et al.ll201 ll) . Both systems also 
display clear evidence for high velocity wings in CO, indicating 
powerful molecular jets. Hence, for these objects, CN may come 
from the outflows rather than from a circum-stellar/binary disk. 
As the line width are however limited, it is also possible that a 
substantial fraction of CN comes from the extended envelope in 
which these objects are embedded. 

Fo r DG Tau, although modeling performed by ISalter et all 
d201 ll) shows the CN line intensity can be produced by a large 
Keplerian disk, we suspect that a similar situation occurs be- 
cause of the strong outflow and complex, clearly non Keplerian, 
circumstellar environment (see Appendix lAl for more details on 
this source). 

For other sources, a disk origin for the CN emission is in- 
ferred from the following qualitative arguments: 
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Fig. 11. CN line width (km.s -1 ) as a function of disk inclination. 
Sources with unknown inclinations are in blue: we used i = 60° 
as the most likely value. Sources with outflows were excluded. 



- The line-width Av derived from the Gaussian fit is of the or- 
der of 2 - 4 km.s -1 , as expected for a Keplerian disk orbiting 
a star of 0.5-2 M Q . 

- Larger line widths are found for the more inclined objects 
(see FigfTTT). This figure does not correct for stellar mass: 
doing so would improve the correlation between line width 
and inclination. 

- The CN outer radius is within 300-600 AU, in good agree- 
ment with the range of interferometrically derived radii in 
previously imaged Keplerian disks. 

- In objects which have already been imaged in CN, MWC 
480, DM Tau, LkCa 15 dChapillon et alj [2012h, GM Aur and 
AA Tau dOberg et al.l2010h . CN unambiguously comes from 
the disk, even in the confused case of AA Tau. 

Note that none of the main CN line profiles appears double 
peaked (e.g. compare CN and 13 CO in MWC 480 Fig|2]l. This is 
a result of the overlap between the three strongest J=5/2-3/2 hy- 
perfine components which are separated by 0.77 and 1 .50 km.s -1 
respectively. These separations are smaller than the overall width 
due to Keplerian shear, significantly blurring the expected dou- 
ble peaked profiles emerging from such disks in (spatially) unre- 
solved observations. However, the velocity separations are larger 
than the local line width (0.3 - 0.5 km.s -1 ), so that in resolved 
(interferometric) images, the patterns from each of the three 
components are spatially separated. Furthermore, the strongest 
J=3/2-l/2 component, which is separated by 17 km.s -1 from the 
nearest one, is only 1.6 times weaker than the J=5/2-3/2 F=5/2- 
3/2 component. In practice, the hyperfine structure of CN is only 
a minor complication in the derivation of the kinematics from in- 
terferometric measurements, especially at the high densities en- 
countered in disks, which ensures that all hyperfine components 
will have the same excitation temperature. 

To the 24 sources detected in CN in our sam ple, we 
sho uld add 5 more disks: GM Aur dOberg et al.l 120101) . LkCa 
15 dChapillon etal.1 1201 2h . CI Tau, CY Tau, and GO Tau 
dGuilloteau et al.ll2013l in prep), leading to a total of 29 detec- 
tions out of 47 objects. 3 (Haro 6-10, T Tau and DG Tau) are 
very unlikely to be disks, and, on the basis of the existence of 
outflows and SO detection, 2 more are perhaps outflows (HL 
Tau and CB26). This still leads to a disk detection rate of at least 
24 / 47, i.e. > 50% at our sensitivity level of ~ 0.1 Jy.km.s -1 . 
Furthermore, in most disks, the CN N=2-l intensity is within a 
factor 2 of that of 13 CO J=2-l (see Figg). 

This demonstrates the ability of the CN 2-1 line to be a good 
tracer of disk kinematics for stars still embedded in molecular 
clouds whose CO and 13 CO line emission could confuse the disk 
emission profile. 



4.2. Molecular disk properties 

Our study is the first large, sensitive survey for 13 CO, CN, ortho- 
H 2 CO, and SO towards disk around Class II PMS stars. Detailed 
sensitive studies of several molecules concerned only a few ob- 
jects: hiJauru^rilyGGTauJ^ 

480 (Dutrev et alJll997t iGuilloteau et al.l[l999t iThi et al.ll20a 



iPietu et all 120071: IChapillon et all 120121: iDutrev et al.l 120111 
Previous multi-molecule surveys were limited to a dozen ob- 
jects, and in general biased to wards the brig htest disks, either in 
CO ( e.g. the DISCS su rvey of lOberg et aLlhoill) , or in contin- 
uum (Salter et al. 2011, used a minimum 1 .3 mm continuum flux 
of 75 mJy). We reach here a sensitivit y level (~ 0.1 Jy. km.s -1 ) 
comparable to or better than that of Oberg et al] d201 lb . and 5 
to 10 times better than that of Sal ter et al.l d201 ll) . Furthermore, 
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Table 5. CN outer radii and line detections 



Source 


Radius 


To 


Line width 


Detection ? 






(AU) 


(K) 


km.s 


i 


CO 


H 2 CO 


SO 


FNTau 


< 170 


15 








Cloud 




CWTau 


< 170 


15 






Disk 






CIDA-1 


< 170 


15 












BPTau 


< 250 


30 












DETau 


< 250 


15 












RYTau 


310 


30 


0.82 ± 


0.27 


Disk 






TTau 


1230 


30 


2.96 ± 


0.02 




Yes 


Yes 


Haro 6-5B 


310 


15 


3.12 ± 


1.17 




Yes 




FTTau 


310 


15 


1.69 ± 


0.34 








DG Tau-B 


450 


15 


2.12 ± 


0.30 




Yes 


Yes 


DGTau 


700 


30 


2.11 ± 


0.09 




Yes 


Yes 


Haro 6-10 


660 


30 


1.50 ± 


0.04 




Yes 


Yes 


IQTau 


560 


15 


4.92 ± 


0.53 


Disk 






LkHa 358 


< 240 


15 






Disk 


Yes 




HH 30 


< 270 


15 






(Disk) 


Yes 




HLTau 


250 


30 


1.79 ± 


0.44 




Yes 


Yes 


HKTau 


320 


30 


8.56 + 


1.96 








Haro 6-13 


< 280 


15 








Yes 




GGTau 


490 


30 


2.28 ± 


0.32 


Disk 


Yes 




UZ Tau E 


310 


30 


6.58 ± 


1.23 


Disk 






04302+2247 


190 


30 


2.99 ± 


0.47 


Disk 


Yes 


Yes 


DL Tau 


560 


15 


1.96 ± 


0.16 




Yes 




DM Tau 


610 


15 


1.44 ± 


0.11 


Disk 


Yes 




AA Tau 


300 


15 


0.95 ± 


0.26 


Disk 






DN Tau 


490 


15 


3.21 ± 


0.42 




Yes 




DO Tau 


310 


15 


0.50 ± 


0.17 


Disk 


Yes 




HV Tau C 


310 


30 


5.97 ± 


0.87 


Disk 






Haro 6-33 


260 


30 


3.05 + 


0.68 




Yes 




DQ Tau 


< 210 


15 












DR Tau 


< 420 


15 












DS Tau 


< 310 


30 








Yes 




UY Aur 


< 300 


30 






Disk 


Yes 




AB Aur 


< 250 


30 






Disk 


Yes 




SU Aur 


< 140 


30 






Disk 






MWC 480 


520 


30 


4.54 ± 


0.34 


Disk 






CB 26 


350 


15 


2.92 ± 


1.00 




Yes 


Yes 


CIDA-8 


<290 


15 








Cloud 




CIDA-11 


<540 


15 












RWAur 


<280 


30 












MWC 758 


<450 


30 






Disk 






CQ Tau 


< 190 


30 






Disk 


Yes 




HD 1632960 


760 


30 


4.16 ± 


0.11 


Disk 


Yes 





Notes. The assumed distance is D=140 pc, except for HD 163296 D=120 pc). For 13 CO, sources in which the emission partially originates from 
a disk are indicated. For H^CO, no distinction is made upon the origin of the emission, apart from the two likely cases of cloud confusion FN Tau 
and CIDA-8. 



the smaller beam size of the 30-m allows to better isolate confu - 
sion cases compared to the JCMT study of Salt er et alJ (1201 ll) . 
The improved sensitivity allows us to detect many more sources, 
such as FT Tau, IQ Tau, V806 Tau (Haro 6-13), UZ Tau, DN Tau, 
DO Tau in the Salter et al. sample. Excluding sources detected in 
SO which are presumably contaminated by outflows, we have a 
total of 32 sources^ (out of 47) which exhibit CN, 13 CO or H 2 CO 
emission from a disk, i.e. 68 % of the sample. 

Figures [12] and Q~3] explore possible correlation of the CN 
intensity with spectral type of the star and 1.3 mm continuum 
flux of the disk respectively. Figure [14] considers the H2CO/CN 
intensity ratio vs the stellar spectral type. Despite the large sam- 



4 This includes the 16 sources detected in CN that have no SO emis- 
sion, plus the Butterfly star, plus 8 sources in 13 CO (including HH30) 
and 2 sources in H 2 CO (Haro 6-13 and DS Tau), to which we add the 5 
sources detected in CN by other studies 



pie size, no correlations are apparent. This lack of correlation is 
also observed for other stellar parameters (mass loss, X-ray lu- 
minosity). However, a number of properties emerge from these 
correlation plots. 

First, CN line is relatively weak (in fact undetected so far) 
in the hottest stars in our sample (SU Aur, CQ Tau, MWC 758, 
AB Aur), although MWC 480 and HD 163296 exhibit strong 
CN emission (FigfT2l). The rather poor molecular content of AB 
A ur has been shown to b e a consequence of its strong UV flux 
bv lSchrever et a l. (2008). The UV flux penetration may further- 
more be enh anced by the larg e inner cavity (radius ~ 100 AU) 
of AB Aur dPietu et al.l l2005h . A similar situation may occu r 
for MWC 758 due to its - 40 AU cavity dlsella et al.l l2010b). 
Chanillon et all (120081 1201 Oh argue that CQ Tau is gas deficient. 
Little is known about the gas disk of SU Aur, but the strong 
X-ray luminosity of SU Aur may play a role in lowering the 
molecular abundances. MWC 480 and HD 163296 have strong 
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Fig. 13. As FigQ~2] but for CN line flux (main fine structure 
group) vs 1.3 mm continuum flux (in mJy). 



CN. Contrary to AB Aur, CQ Tau and MWC 758 which are Type 
I sources , these two sources a re Type II according to the classifi- 
cation of iMeeus et al.l (1200 ll) . which can indicate self-shadowed 
disks. Such a configuration would seriously limit the UV pene- 
tration, leading to a chemical composition much more similar to 
those of T Tauri disks. 

Second, at the level of our sensitivity, we do not detect the 
CN line for spectral types cooler than Ml. 5 (Figfl2l. Assuming 
similar CN column densities than in the warmer stars, our sensi- 
tivity translates into upper limits of order 150 - 250 AU for the 
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Fig. 14. Correlation between the H 2 CO to CN flux ratio and stel- 
lar spectral type. Red indicates upper limits (no H2CO), black 
lower limits (no CN). H 2 CO in FN Tau and CIDA 8 is likely to 
be cloud contamination. Outflow sources are in italics. 



disk outer radius. Thus, our negative result may only be the con- 
sequence of smaller disk sizes for lower mass stars, rather than 
implying a change in chemistry at this spectral type threshold. 

Third, although CN intensity does not correlate with the con- 
tinuum flux at 1.3 mm (Fig[T3l. which is a tracer of the dust 
mass, there is a clear threshold. Only 4 sources with S v (1.3 mm) 
> 60mJy are undetected in CN (CW Tau, CQ Tau, DQ Tau, and 
DR Tau which has been observed with relatively limited sensi- 
tivity) while only 4 sources are detected among the 17 sources 
with S v ( 1.3 mm) < 60mJy. 

Comparing the emissions from CN and H2CO, we find that 
to first order, the H2CO / CN ratio divides the sources in two cat- 
egories (see FigfT4b. [H2CO/CN]> 0.5 is found in warm sources 
(Herbig Ae stars) or outflow/envelope dominated objects, while 
it is below 0.3 for the disk sources. This trend is consistent 
with our current understanding of chemist ry. In disks, H 2 CO 
is for med in the gas phase and on grains (iSemenov & Wieba 
l20lTh . Therefore, H 2 CO can be released in the gas phase at high 
enough temperatures, either because of large stellar flux (HAe 
stars) or because of shock heating (outflow sources). Enhanced 
abundance of H 2 CO has been observed i n many outflows (e.g. 
iBachiller et ail 120011: iTafalla et alj l2010h . As a consequence, 
sources with weak or no H2CO emission mostly correspond to 
cold objects without outflows. While CN appears as a clear disk 
tracer in the M1-K5 spectral type range, the situation is more 
complex for the HAe stars, for which H2CO is so metimes a bet- 
ter pro be than CN, but remains weak. Using the Baraffe et al. 
( 1 998) tracks, the M 1 limit corresponds to about 0.5 M Q for stars 
in the 1 to 3 Myr age. 

Sulfur-bearing mo lecules (CS, SO, H?S and CCS) have been 
recent ly searched by iDutrev et al.1 d201 lb and fChapillo n et al.l 
(120121) in the four disks surrounding MWC 480, LkCa 15, GO 
Tau and DM Tau. So far, only CS has been detected, neither SO 
nor H2S. Our survey provides a much larger sample for SO, and 
confirms its general weakness in disks (only one source detected, 
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Fig. 15. Histogram of detections as a function of 1.3 mm con- 
tinuum flux range. Filled: disks, hatched: outflow contaminated 
sources, empty: non detections. 



04302+224 2, our results cont radicting the detection in AB Aur 
reported bv lFuente et al.1 (120101) '). On the other hand, we find SO 
to be ubiquitous in outflow sources. The SO chemistry in disks 
is not fully unders tood: current model s predict too much SO 
and H2S molecules iDutrev et al.l (1201 ll) . most likely as a result 
of incomplete treatment of Sulfur chemistry on grain surfaces. 
However, SO is a clear tracer of sh ocks, with high ab undances 
in extremely high velocity outflows (Tafall a et al.l20ldl) . 

The absence of correlations in Fig fT2lfT3l justifies a posteri- 
ori our simple source-independent hypothesis (essentially con- 
stant CN surface density) to interpret the line emission. In this 
framework, the CN flux is a simple tracer of disk size, and our 
study shows that at least 50% of the young stars are surrounded 
by large (> 300 AU) disks, in our specific case of isolated star 
formation. The overall disk detection rate is displayed in Fig[T5l 
Hatched areas are for sources with molecular outflows (traced by 
SO), filled areas are disks detected in CN, 13 CO or H 2 CO, and 
empty histograms non detections. The distribution of disk sizes 
exhibits a strong peak near 300 AU; however, because of sen- 
sitivity limitations, our survey becomes incomplete for smaller 
disk sizes. 

It is beyond the scope of our angularly unresolved data to 
provide detailed modeling of the disk emission for individual ob- 
jects. The above considerations on the detected lines are broadly 
consistent with prior knowledge of disk chemistry. The large CN 
detection rate supports the simple idea that CN is present only in 
a photodissociation layer above the disk plane, with a relatively 
uniform surface density over a large range o f source properties 
(|JonkheidetalJl2007t IChapillon et al1l2012h . The lower detec- 
tion rate for sources with low 1.3 mm flux may be the result of 
a change of chemical regime for disks with low mass. However, 
it is also consistent with a simpler idea, namely that these faint 
disks are just much smaller than the others. Higher sensitivity, 
spatially resolved observations can solve this issue. 



and C 17 rotational lines near 206 - 228 GHz with the IRAM 
30-m telescope. It reveals that : 

- 13 CO is strongly affected by confusion, either from molecu- 
lar clouds, or from outflows and envelopes. 

- CN is a good, nearly confusion free, tracer of disks for 
stars in the M1-K5 spectral type range (i.e. stars more mas- 
sive than about 0.5 M according to the Bar affe et alj d!998h 
tracks). 

- The warmest stars exhibit two very different behaviours: 
some have strong CN, but most have undetectable CN and 
weak H2CO emission. 

- SO is ubiquitously found in outflow-driving, embedded 
sources, but exceptional in disks (only 1 source). When SO 
is present, H2CO is stronger than CN. 

- CN N=2-l is strong enough to be a good tracer of disk kine- 
matics in at least half of the objects. 

- The disk detection rate in CN lines exceeds 50 %. The de- 
tected disks are large, > 300 AU in radius. 

- The overall gas disk detection rate (by any molecular tracer) 
is a: 68 %. However, gas disk detection is much more dif- 
ficult for sources with small (< 60 mJy) 1.3 mm continuum 
flux. We cannot distinguish at present whether this is due to 
smaller disk masses or smaller disk sizes (< 200 AU radius). 

Our results only apply to regions of isolated star formation. In 
more crowded regions (e.g. p Oph), confusion will unavoidably 
be larger. The disk size distribution may also be different. It 
is however likely that CN will remain a better tracer of disk 
kinematics than the much more confused I2 CO or 13 CO, or the 
significantly weaker H2CO emission. Looking to the advent of 
ALMA, with 8 times greater collecting area than the IRAM 30- 
m, better antenna efficiency, and expected system temperatures 
around 100 K at this frequency, will provide a sensitivity gain 
of order 16. Thus ALMA should be able to obtain, in the same 
integration time, a 4 times better limit on the outer radius. This 
may open for the first time the possibility to probe beyond the 
Ml limit, unless the chemistry changes to reduce drastically the 
CN abundances around the lower mass stars. 
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5. Summary 

We have performed a sensitive survey of 42 young stars probing 
isolated star formation regions in CN, ortho-H2CO, SO, 13 CO 
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Appendix A: Comments on individual sources 

FN Tau: the tentative detection of H2CO is at an unusual veloc- 
ity of -0.5 km.s -1 . 



reliable. Our CN line width, though unprecise, suggests a higher 
stellar mass. The asymmetric H2CO spectrum may be affected 
by confusion in the blueshifted part. 



CW Tau is dominated by confusion (see discussion in Appendix 
ICt . Despite strong confusion, line wings of 13 CO could in- 
dicate a disk at velocity 6.2km.s~ 1 and with a line width of 
3.3 + 0.4 km. s -1 . CW Tau drives a bipolar outflow that can be 
traced in o ptical forbidden li nes over 4" - 6" along both outflow 
directions (Iffirth et al.ll 1994 . 



CI DM lBricenoetalJ(ll993h is a late spectral type object (M5.5 



White & Basri 



2003), with weak 1.3 mm continuum emission 



(ISchaefer et alJl2009h . The line of sight towards CIDA 1 inter- 
cepts the same molecular cloud than that of CW Tau and is 
highly confused. The tentative fit of the 13 CO line wing is not 
conclusive on the presence of a disk. 



BP T au harbors a small (120 AU radius) disk dDutrev et alJ 
12001 and strong confusion in 12 CO. No signal at all is detected 
in these new observations, but the sensitivity towards this source 
is lower than average. 



DETau has relatively weak 1.3 mm continuum (36 mJy). There 
is some confusion in 13 CO, but no evidence for disk emission. 



RY Tau: the CN detection is marginal (3.5 <x). A fit of a 3 km.s -1 
wide gaussian to the 13 CO spectrum also indicate a potential 
disk, at the 5cr level, but with a different veloci ty. RY Tau is 
a susp ected close binary, with sep aration > 3 AU (Bertout et alJ 



1999). Velocity measurements by Ngu yen et alJ d2012 ) indicate 
a possible SB1 with radial heliocentric velocity in the range 
16.7 - 19.4km/s, which is equivalent to 7.7 - 10.4km/s in the 
LSR, and more consistent with the 13 CO velo city. High reso- 
lution image at 1.3 mm bv llsella et all i20l0eb indicate an in- 
clination of 66° and a morphology consistent with a deficit of 
emission in the inner 15 AU. 



T Tau is a triple system, with a compact (~ 60 AU) c i rcum- 
stellar disk around the sing le T Tau N (Ake son et alJ Fl998: 
iGuilloteau et al.ll201 lh . Its molecular environment is complex, 
and all the detected lines are most likely coming from an outflow 
or heated e nvelope, perhaps the o utflowing cavity walls imaged 
in 13 CO bv lMomose et all dl996h . 



Haro 6-5 B, also known as FS Tau B, is a nearly edge-on, 
deeply embedded object. Confusion in 13 CO is very strong. 
Previous interferom etric measurements of 13 CO J=l-0 revealed 
a complex region dDutrev et all [1991 . HST imag es indicate 
a 300 AU radius disk at high inclination, 70-80° dKrist et alJ 
119981: |Padgett et al.llT999h . FS Tau B possess a jet an d a coun- 
ter jet , as well as an optically visible cavity wall dLiu et alJ 
120121) . Our combined detection of CN and H2CO may indicate 
a compact, warm disk or some emission from an outflow or 
envelope. 5" resolution images in 13 CO J=l-0 with the NMA 
dYokogawa et all l2002h revealed a disk at V LSR ~ 7.0km.s _1 , 
with a kinematics dominated by rotation. They indicate motions 
compatible with Keplerian motions around a central star of 0.25 
M only. Given the high confusion level, this value may be un- 



FT Tau is a deeply embedded, low luminosity star with un- 
kn own spectral type . A c ompact (60 AU) dust disk was resolved 
by IGuilloteau et alJ d201 lh . who also detected 12 CO emission, 
suggesting a small disk at low inclination around a 0.7 -1.0 M Q 
star. The clear detection of CN confirms the disk interpretation. 
13 CO emission from the disk may be present at redshifted veloc- 
ities compatible with CN, but confusion makes any quantitative 
estimate impossible. 



DG Tau is a well studied object, driving an optical microjet. The 
sys tem inclination remain unclear: between 38 and 45° from the 
jet dPvo et al.l2003fclEisloffel & Mundtl 19981). and 20 to 32° from 
the dust disk dlsella et al|201 Oat IGuilloteau et alj|201 lh . Images 
in 13 CO (2-1) dTesti et alJ2002l) indicate a complex morphology, 
a large fractio n of the emission co ming from a remnant enve- 
lope (see also iKitamura et alii 19961) . There is no clear evidence 
for Keplerian rotation. The relatively narrow lines of H2CO and 
SO could originate fro m this remnant en velope. CN(2-1) was 
previously detected by Salte r et al.1 (Hoilfl, who modeled this 
as coming from a Keplerian d isk. Our CN line w idth is consis- 
tent with the measurement of ISalter et al. (2011), but the large 
derived outer radius, combined with the complex morphology, 
suggests CN is most likely coming from the outflow or the rem- 
nant envelope. 



DG Tau B is a totally o bscured object , driv ing a powerful one- 
sided molecula r outflow [ Mitchell et al. ( 199 7|). High inclination 
(64 - 75°, see IGuilloteau et al.1 1201 It lEisloffel & Mundtl fl998) 
however yield moderate line-of-sight velocities for the outflow. 
The rich spectrum may originate from the outflow, but CN lacks 
the narrow component detected in H2CO and C n O, which is 
presumably due to the molecular cloud. 



Haro 6-10 (also known as GV Tau) is a close binary, with 
two compact (< 30 A U), optically thick, circumstellar disks 
dGuill oteau et aT] |201 lh . but no known circumbinary disk. The 
circumstellar disks are misali gned and the system is e mbed- 
ded in a common envelope (Roccatagliat a et alJ I201 lh . H2O 
maser emission has been recently reported, an d suggest that 
the S outhern component is itself a close binary (Wilkin get al.l 
l20ll . As for DG Tau B, the CN line is much wider than those 
of H2CO and C n O The preeminent red wing in 13 CO is likely 
to trace the outflow. 



IQ Tau was observed in 12 CO by ISchaefer et ail ((2009), who 
did not detect any significant emission with the IRAM Plateau 
de Bure interferometer at 5" resolution. The CN emission ex- 
hibits the clear signature of Keplerian rotation (although the ex- 
pected double peaked aspect is reduced because of the hyperfine 
component blending). CN is apparently stronger than 13 CO in 
this source. A tentative fit of a disk component to the 13 CO spec- 
trum suggest the redshifted emission may still be masked by the 
cloud. 



Although their Figs. 10 and 12 mark it as the 3-2 transition 
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LkHa 358: I3 CO is heavily confused. A small molecu- 
lar disk was repor ted from interferometric measurements by 
Sch aefer et al.1 (120091) from 12 CO emission. An attempt to fit a 
disk component to the 13 CO line wings is consistent with this 
interpretation, provided 13 CO has significant optical depth (> 2 
on average). The narrow emission from C I7 most likely comes 
from the surrounding cloud, as the ~ 3.5<x signal from H2CO. 

HH 30: iPetv et al.1 d2006l) had previously imaged the 13 CO 
emission from this emblematic edge-on object. The narrow line 
from C n O likely comes from the cloud, as the ~ 3cr signal from 
H 2 CO. 

HL Tau is the archetype of the deeply embedded, outflow driv- 
ing sourc es. The redshifted side of the outflow is known to be 
brighter (ICabrit et al.l Il996t iMoninet ail 1 19961) . From the line 
profiles, C I7 0, H2CO and especially SO (which is significantly 
redshifted), must come from the outflow. However, the CN line 
is much narrower and exhibit a systemic velocity of 4.7 km.s -1 , 
and thus could be a pure disk tracer. 

H K Tau is a binary star, with o ne component exhibiting an edge- 
on disk dMcCabe etal.ll2011l) . The relatively broad CN line is 
consistent with emission from this edge-on disk around HK Tau 
B. The blue shifted side of the disk may also have been detected 
in 13 CO. 

Haro 6-13, also kno wn as V 806 Tau, was imaged in 12 CO by 
Scha efer et al.l (I2009I) . who based on agreement in position angle 
between CO and dust emission, attributed the emission to a disk 
of at least 180 AU rad ius. HCO + 3-2 emission has been detected 
by ISalter et all (1201 ll) . The H2CO spectra are complex, as the 
C 17 one, and perhaps result of superposition of relatively broad 
emission with a narrow (0.5 km.s -1 ) absorption feature. 

GG Tau was included as a reference sourc e in our sample. Ou r 
results agree with the initial detections of iDutrev et al.l (Il997l) . 
H2CO is well detected in the combined spectrum. 

UZ Tau E is a spectroscopic binary part of the UZ T au hierarchi- 
cal qu adruple system. The CO disk was detected bv lSimon et al.l 
(2000), who found an inclination of 55° and an outer radius of 
300 AU. Our derived radius from CN is in excellent agreement 
with this value. 13 CO emission from the disk is also detected, in 
the line wings, but still exhibits confusion. 

IRAS04302+2247 (the Butterfly star iWolf et all l200l . is an 
edge-on system. Recent IRAM interferometric measurements 
indic ate a stellar mass o f order 1 .6 M and a disk radius of 260 
AU (IDutrev et al.ll2013"l in prep). Although this source exhibits 
a large H2CO / CN ratio, which we in general attribute to out- 
flows, the H2CO and SO emissions come from the disk. From 
the hyperfine ratios, CN appears potentially optically thick in 
this source, in which case the emission must be dominated by 
smaller radii. The rather large linewidth concurs with this inter- 
pretation. 

DL Tau was imaged in 12 CO by ISimon et all (120001) who 
pointed out the strong confusion making the inclination deriva- 



tion unreliable. Confusion is also large in 13 CO, masking all the 
emission except for the red-shifted wing. Like IQ Tau, DL Tau 
shows strong, double peaked line profile in CN. The outer ra- 
dius derived from CN is consistent with the value quoted by 
ISimon et all (l2000h for CO. H 2 CO is also detected. 

A A Tau: CN is clearly detected, and weak emission is also vis- 
ible in 13 CO, despite confusion. However, the disk con tribution 
in 13 C O is much smaller than the intensity reported bv lGreavesI 
d2005l) (0.6 K.km.s 1 at JCMT, i.e. 12 Jykm.s" 1 ), which must 
have included clou d emission. CN in AA Tau was discovered by 
IQberg et aD (l2010l) . We recover only about half of their quoted 
flux (1.7 Jykm.s -1 ), and our derived velocity suggests that the 
blue-shifted part of the CN spectrum is masked by confusion, 
like the 13 CO one. 

DN Tau is very similar to DL Ta u and IQ Tau, with CN stronger 
than 13 CO. ISchaefer et al.l d2009l) failed to detect 12 CO from this 
object with the IRAM interferometer, again presumably because 
of strong confusion with the cloud. 

DO Tau shows clear CN detection (4<x), but the line width is 
poorly constrained. The line can be narrow, so (positive) con- 
tamination by the cloud is not fully excluded. Identification of 
a disk component in I3 CO is difficult, as the redshifted side is 
dominated by the cloud. 

HV Tau C, like HK Tau B, is a compact (80 AU radius) edge on 
disk in a multiple system, here a hie rarchical triple. J^CO emis- 
sion from the disk was resolved by iDuchene et al.l yOlO) with 
the IRAM array. The broad CN line (AV = 6km.s _1 ) is con- 
sistent with an origin in a compact disk. The CN flux suggests 
a somewhat larger radius (300 AU in Table [3), or equivalently, 
non negligible optical depth. The observed hyperfine ratios are 
indeed compatible with substantial optical depth, r = 7.0 + 2.7, 
which would also yield a lower line width, 4.2 + 0.6km.s~'. 
I3 CO emission from the disk produces the line wings visible in 
FigE26] 

Haro 6- 33 (also known as IR AS04385+2550) was detected in 
I2 CO by Schaefer et al. (2009), who interpreted the emission as 
coming from a Keplerian disk, despite strong confusion. There is 
also substantial cloud contamination in I3 CO, but CN and H2CO 
are detected. The derived disk radius in CN is in good agreement 
with the 12 CO result. 

DQ Tau exhibits huge confusion, and no detection at all. 

PR Tau was rep orted as tentatively detected in CN(2-1) by 
ISalter et ail d201 lb at a velocity of 1 1 km.s~' with an integrated 
flux of 3.5 ± 1.3 Jy km.s~'. Our spectrum shows no line at all, 
with a 3<r upper limit on the flux of 1.4 Jy km.s -1 . Confusion 
with the cloud is huge in 13 CO: only the red wing between 10 
and 1 1 km.s~' may be attributable to a disk component. 

DS Tau shows no apparent confusion. A weak (4cr) line is de- 
tected when combining both 0-H2CO transitions, with a veloc- 
ity (6.0 km.s -1 ) and linewidth (3 km.s -1 ) consistent with disk 
emission, but no CN or I3 CO. 
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UY Aur is a 0.8" separation binary system. The 13 CO line is 
strong, and consistent w ith an origin i n the complex circumbi- 
nary structure imaged bv lDuvert et al.l d2000t) . There is a Act de- 
tection of H2CO with a compatible kinematics. 

AB Aur: the strong I3 CO line ema nates from a co mbination 
of the disk/ring structure imaged by iPietu et a"fl (120051) and the 
surrounding envelope. The clear detection of C I7 confirms 
the large CO content of this sour ce. H2CO is convin cingly de- 
tected, confirming the result of iFuente et al.l (|2010) using the 
para-H^CO 3o,3 - 2o,2 line at 218.2 GHz. However, we find no 
evidence at all for CN or SO from the disk, which were claimed 
to be detected by these authors. We only have a narrow (0.4 
km.s -1 ), 3.5cr feature in CN(2-1) at a velo city consistent with 
cloud emission. The CN ( 1 -0) line quoted bv lFuente et al.l {20 10) 
must have been residual emission from the cloud. For SO, they 
report a line intensity of 26 mK km.s~' for the 34 - 23 line, and 
60 mK km.s -1 for the 56 - 4g transitions, while we obtain a 3cr 
upper limit of 21 mKkm.s 1 for th e 54 - 4^ using the A B Aur 
velocity and line widths. Note that IFuente et al.l d2010l) used a 
larger throw (120") than us in the wobbler switching, so that 
their results are more susceptible to residual emission from the 
envelope. 

SU Aur: despite confusion with the cloud, the 13 CO spectrum 
exhibits a broad component compatible with emission from 
a disk. No ot her molecule is detect ed. SU Aur is a strong 
X ray source iFranciosini et al.l (|2007). Near IR imaging indi- 
cates a highly inclined disk (^ 65°) and ou tflow cavity walls 
(lAkeson et : alJl2002t IChakrabortv & Gdl2004h . 

MWC 480 was included here as a reference source in our sam- 
ple. It has been e x tensively studied with the IRAM mm array 
IPietu et al.1 (120071) : IChapillon etail (|2012|) . We find the double 
peaked 13 CO profile characteristic of an isolated Keplerian disk. 
CN emis sion is in excellent agre ement with the interferometric 
results of IChapillon et al.l (120121) . The detection of C 17 for the 
first time in this source confirms the significant CO abundance in 
this object. Note that, as expected, the C n O profile is apparently 
larger than the 13 CO one. 

CB26 is a deeply embedded, edge-on object driving an unusual 
outflow Launhardt et al .1 (120091) in an isolated Bok globule. Like 
in HL Tau, we detect here strong H2CO and SO emission, but 
only a weak CN line. If interpreted by disk emission, the CN 
line is consistent with a disk size about 250 AU. 

CIDA 8 has weak continuum e mission and no CO disk from 
interferometric measurements of Sch aefer et al.l d2009l) . Cloud 
contamination exists, but appears limited to a very narrow ve- 
locity range. There is a surprising detection of H2CO just outside 
the range visibly contaminated in 13 CO, but with a narrow line 
width which suggests it originates from a cloud. 

CIDA 1 1 shows no d etection at all. I t was also not detected in 
continuum and CO bv lSchaefer et al. (2009). 

RW Aur is a bin ary system that was imaged in 12 CO by 
ICabrit et"aTI ([2006), who find evidence for a very small (50 AU 



radius) disk, but did not detect 13 CO emission. Our negative re- 
sult is consistent with their findings. 

MWC 758 is a 1.8 M star surrounded by a di sk where CO is 
appar ently depleted and mostly optically thin (Cha pillon et al.l 
2008), despite high temperatures (30 K). We find a relatively 
weak line of 13 CO, consistent with this interpretation. No other 
molecule is detected: the marginal signal in C n O is not at the 
expected disk velocity. The li mited CO opacity c ould be related 
to the large cavity detected bv llsella et al.ld2010bl) in continuum. 

CQ Tau is a "twin" of MWC 758, also exh ibiting a warm 
(> 50 K) disk with optically thin CO emission (IChapillon et al.l 
2008) , but its continuum emission is strong ly centrally peaked 
dBanzatti et al.ll201 it [Guilloteau et al.ll201 II) . Although different 
distances are often quoted based on Hipp arcos par a llax, these 
differences are not significant and IChapillon et al.l (120081) ar- 
gue that both CQ Tau and MWC 758 are also located at 140 
pc. We also detect 13 CO with a similar intensity than in MWC 
758. H2CO is detected at the 5<x level, with kinematic param- 
eters quite consistent with those derived from 13 CO. However, 
we find no evidence for CN. Our 3 cr limit of 0.28 Jykm.s 1 
is marginal ly consistent with t he flux of 0.22 ± 0.11 Jykm.s 1 
reported bv lOberg etail (l2010l) . 

HP 163296 is an isolated Herbig Ae star located at 120 pc 
(Ivan Leeuwenl 120071) . I t harbours a large (~ 550 AU) CO disk 
inclined at about 45 ° (llsella et all 12007). and has been exten- 
sively studied in CO isotopologu es, includ i ng C 1 7 Q dOi et al.l 
1201 ll) . CN N=3-2 was detected bv lThi et al.l d2004l) . The strong 
CN N=2-l emission is compatible with an optically thin line. 
The outer radius in Table [5] is slightly larger than that mea- 
sured from CO. The discrepancy may have several origins. First, 
as the st ar is very luminous, CN may be in a warmer region. 
Second, Hugh es et al.l d201 ll) report a rather large turbulence 
level (^ 0.3km.s~') from CO measurements, larger than as- 
sumed. We also report the first detection of H2CO, through 
the ortho-H2CO 3n - 2i2 transition. C omparison of o ur 13 CO 
spectrum to the interferometric result of lOi et all (1201 ll) suggest 
some contaminatio n near 2 - 4km .s -1 (and also 13km.s _1 as in 
the 12 CO spectra of lOi et al.1 (120111) 1. 

Appendix B: Spectra for individual sources 

This appendix displays the spectra towards the various sources. 
For each source, the top panel shows the 13 CO J=2-l spectrum. 
The continuous red line, if present, is a Gaussian fit obtained af- 
ter masking the confused velocity range. The bottom panel dis- 
plays on a common scale from top to bottom the spectra of: I3 CO 
J=2-l (with fit as in top panel), C n O J=2-l, CN N=2-l, ortho 
H2CO 3 13 - 2i2, ortho H2CO 3i2 - 2n, the average of both o- 
H2CO transitions, and SO 54 - 43, arbitrarily shifted in intensity 
to avoid overlap. The intensity scale is antenna temperature (T^ 
in K): conversion to flux density can be obtained using a factor 
of9Jy/K. 
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Fig. B.l. Spectra of the observed transitions towards FN Tau Fig. B.3. Spectra of the observed transitions towards CIDA-1 
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Fig. B.2. Spectra of the observed transitions towards CW Tau Fig. B.4. Spectra of the observed transitions towards BP Tau 
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Fig. B.5. Spectra of the observed transitions towards DE Tau Fig. B.7. Spectra of the observed transitions towards T Tau 
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Fig. B.6. Spectra of the observed transitions towards RY Tau Fig. B.8. Spectra of the observed transitions towards Haro 6-5 B 
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Fig. B.9. Spectra of the observed transitions towards FT Tau Fig. B.ll. Spectra of the observed transitions towards DG Tau 
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Fig. B.10. Spectra of the observed transitions towards DG Tau B Fig. B.12. Spectra of the observed transitions towards Haro 6-10 
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Fig. B.13. Spectra of the observed transitions towards IQ Tau Fig. B.15. Spectra of the observed transitions towards HH 30 
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Fig. B.14. Spectra of the observed transitions towards LkHa 358 Fig. B.16. Spectra of the observed transitions towards HL Tau 
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Fig. B.17. Spectra of the observed transitions towards HK Tau Fig. B.19. Spectra of the observed transitions towards GG Tau 
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Fig. B.18. Spectra of the observed transitions towards Haro 6-13 Fig. B.20. Spectra of the observed transitions towards UZ Tau E 
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Fig.B.21. Spectra towards 04302+2247 



Fig. B.23. Spectra of the observed transitions towards AA Tau 
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Fig. B.22. Spectra of the observed transitions towards DL Tau Fig. B.24. Spectra of the observed transitions towards DN Tau 
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Fig. B.25. Spectra of the observed transitions towards DO Tau Fig. B.27. Spectra of the observed transitions towards Haro 6-33 
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Fig. B.26. Spectra of the observed transitions towards HV Tau C Fig. B.28. Spectra of the observed transitions towards DQ Tau 
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Fig. B.29. Spectra of the observed transitions towards DR Tau 
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Fig. B.31. Spectra of the observed transitions towards UY Aur 
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Fig. B.30. Spectra of the observed transitions towards DS Tau Fig. B.32. Spectra of the observed transitions towards AB Aur 
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Fig. B.33. Spectra of the observed transitions towards SU Aur Fig. B.35. Spectra of the observed transitions towards CB 26 
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Fig. B.34. Spectra of the observed transitions towards MWC 480 Fig. B.36. Spectra of the observed transitions towards CIDA-8 
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Fig. B.37. Spectra of the observed transitions towards CIDA-1 1 Fig. B.39. Spectra of the observed transitions towards MWC 758 
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Fig. B.38. Spectra of the observed transitions towards RW Aur Fig. B.40. Spectra of the observed transitions towards CQ Tau 
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Fig. B.41. Spectra of the observed transitions towards DM Tau 



Appendix C: The CW Tau / CIDA-1 line of sight 

CW Tau and CIDA-1, which are only separated by 1 arcmin, ex- 
hibit strong confusion in 13 CO, and are the only sources in which 
confusion in CN is obvious. Conf usion is also visib le in H2CO 
and even C n O. In addition, in the lSalter etal] (1201 lh study, CW 
Tau is the source which has extended HCO + (3-2) emission. 

In the beam-switched observations, the detection of hyper- 
fine components of CN (Fig JC. 1 b suggests significant optical 
depth. Optically thick CN from clouds could introduce difficul- 
ties in using CN as a kinematical disk tracer, as it could absorb 
the disk emission. This potential problem is mitigated by sev- 
eral independent arguments. First, the cloud line width is small, 
so only a small fraction of disk kinematics will be affected. 
Second, CN has several hyperfine components, of very differ- 
ent opacities: this would help in modeling the absorption layer. 
Furthermore, anomalously high (r » 1) apparent opacities can 
be the result of opacity gradients in a cloud with more moder- 
ate (t > 1) average opacity, as the strongest lines saturate first. 
However, observations of the CW Tau line of sight in frequency 
switching shows a flat-topped, 2 K line of 13 CO, with a line 
width of 1.40 km.s -1 and the same velocity than the CN signal, 
but no significant CN emission at a level of 50 mK (3(f). Thus, 
in the beam-switched observations, the detected signal probably 
originates from emission in one or both of the reference beams, 
and in the cloud, CN is most likely not much excited. Although 
high opacities cannot be excluded, the unusual hyperfine line 
ratios may also be related to weak excitation-dependent anoma- 
lies, since CN collision rates display signific ant dependencies 
on the hyperfine levels (Kaluginaet aQ 2012b . Non negligible 
opacities, low excitation and possible hyperfine anomalies have 
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Fig. B.42. Spectra of the observed transitions towards HD 
163296. The observed SO transition is not the same as in all 
other sources. 
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Fig. C.l. Spectra of CN towards CW Tau. The best fit with the 
hyperfine components is drawn below. The upper profile indi- 
cates the expected profile for an optically thin line. Detection of 
even weak components suggests high optical depth (~ 20 - 30). 
The intensity scale is T* A (K). 



alread y been reported for CN in dark clouds by ICrutcher et alj 
(11984 . 



Appendix D: Formation of Lines in Keplerian Disks 

Line formation in a Keplerian disk is stron gly constrained by 
the velocity gradient dHorne & Marshlll986h . We use here this 
property to explore to what extent the detected CN lines trace 
the (putative) circumstellar disks. Unless the disk is seen face 
on, only a fraction of the disk projects at any given velocity. This 
fraction can be estimated wit h simple reasoni n g: we follow here 
the derivation presented by iGuilloteau et al.l (12006). r, being 
the cylindrical coordinates in the disk plane, the line of sight 
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Fig. D.l. Regions of the disk which yield equal projected ve- 
locities v bs- The ellipse is the projection of the disk outer edge 
observed with an inclination of 45°. 



velocity (in the system rest frame) is 



inclination and line width disappears for optically thin lines: the 
integrated line flux scales as the number of molecules. 

From Eq |D.6l we can derive the outer radius of the molecule 
distribution provided with have an estimate of p 



R ml , — D 



§ Sydv 



B v (T )(pAV)ncos{i) ) 



1/2 



(D.7) 



Equation lD.7l is valid for single lines. For molecules with hyper- 
fine structure like CN, line blending must be accounted for. We 
do so in taking into account the relative weight of the hyperfine 
components: the integrated flux of the J=5/2-3/2 fine structure 
group presented in Table [3] is multiplied by 0.45 to obtain the 
flux in the strongest hyperfine component, which then used in 
Eq lD.7l to derive the outer radius. The fitted opacity, if available, 
is used to provide an estimate of p (p — t, t being the sum of the 
opacities of all hyperfine components), which is otherwise set to 
0.3. 



V ohs (r,0)= y/GMJr sin; cos 9 



The locii of isovelocity are given by 
r{6) = (GM,/V 2 bs ) sin 2 /cos 2 9 



(D.l) 



(D.2) 



With a finite local line width Av (assuming rectangular line 
shape for simplification), the line at a given velocity V b s origi- 
nates from a region included between and r s {9) : 



n(9) 

r s (9) = min 



GM t 



(V obs + Av/2) 2 



■ sin 2 i cos 2 9 



GM* 



(V obs - Av/2) 2 



sin 2 / cos 2 9 



(D.3) 
(D.4) 



Figure lDTI indicates the regions of equal projected velocities for 
6 different values: V D b s > v^, V b s = vj, V D b s < vj, and their 
symmetric counterpart at negative velocities, where vj : 



v,/ 



= y/GM t /R 0Ut sini 



(D.5) 



is the projected velocity at the outer disk radius R oa . Figure ID~T1 
shows that the fraction of the disk covered by the gas at velocities 
Vobs — is of order Av/vj, and drops very rapidly for larger 
velocities. The larger area covered for V b s - vj explains the 
classical double peaked line profiles emerging from Keplerian 
disks. 

For moderately inclined (i < 70°) disks in Keplerian rotation, 
the integrated line flux is given by 



Sydv 



B v (T )(pAV)nR 2 o JD 2 cos(i) 



(D.6) 



where To is the disk temperature (assumed uniform for sim- 
plicity), AV the local linewidth, R out the disk outer radius, D 
is the source distance and i the disk inclination, p is a factor of 
order of the line of sight r; for optically thin lines, saturating 
as oc log(T/) for large optical depths (see iGuilloteau & Dutrevl 
1998, their Fig. 4). The optically thick limit is obtained because 
of the Keplerian shears which, as demonstrated above, limits 
the maximum fraction of the disk covered at any velocity to 
6V/vd- 6V is the local linewidth including opacity broadening, 
so 5V ~ log(r/)AV for large opacities. 

Note that as r; t p / cos(/), where t p is the line opacity 
perpendicular to the disk plane, and since t p oc E/AV, where 
£ is the molecule surface density, the (apparent) dependence on 
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